Comparative Bending Dynamics in DNA with and without A-tracts 



Alexey K. Mazur and Dimitri E. Kamashev 
Laboratoire de Biochimie Theorique, CNRS UPR9080 
Institut de Biologie Physico-Chimique 
13, rue Pierre et Marie Curie, Paris, 75005, France. 
FAX:+SS[0]1. 58.41. 50.26. Email: alexey@ibpc.fr 



o 
o 

(N 
Oh- 

o 



Or 
I ' 

o 



C/3 

O 
Oh: 



> 

OO 

m 
o 

OO 
O 



c*3 
O 

Oh. 



X 



The macroscopic curvature of double helical DNA in- 
duced by regularly repeated adenine tracts is well-known but 
still puzzling. Its physical origin remains controversial even 
though it is perhaps the best-documented sequence modula- 
tion of DNA structure. We report here the results of compara- 
tive theoretical and experimental studies of bending dynamics 
in 35-mer DNA fragments. This length appears large enough 
for the curvature to be distinguished by gel electrophoresis. 
Two DNA fragments, with identical base pair composition, 
but different sequences are compared. In the first one, a single 
A-tract motif was four times repeated in phase with the heli- 
cal screw whereas the second sequence was "random". Both 
calculations and experiments indicate that the A-tract DNA 
is distinguished by the large static curvature and character- 
istic bending dynamics, suggesting that the computed effect 
corresponds to the experimental phenomenon. The results 
poorly agree with earlier views that attributed a decisive role 
in DNA bending to sequence specific base pair stacking or 
binding of solvent counterions, but lend additional support to 
the hypothesis of a compressed frustrated state of the back- 
bone as the principal physical cause of the static curvature. 
We discuss the possible ways of experimental verification of 
this hypothesis. 



INTRODUCTION 

It is generally accepted that the base pair sequence 
can affect the overall form of the DNA double helix. In- 
trinsic DNA bending is the simplest such effect. Natural 
static curvature was discovered nearly twenty years ago 
in DNA containing regulas repeats of A n T m , with n + 
m > 3, called A-tractsatali Since then this intrigu- 
ing phenomenon has been actively studied, with sev- 
eral rffafrirjrid reviews of the results published in different 
yearsU'ErQ'Q'ErO. It is known that the curvature is directed 
towards the minor grooves of A-tracts and/or the major 
grooves of the junction zones between them, and that its 
magnitude is around 18° per A-tract. However, the exact 
sites and the character of local bends remain a matter of 
debate as well as their mechanism and physical origin. 

Already the pioneering conformational calculations of 
the seventies showed that the DNA double helix exhibits 
significant, bf usability which is anisotropic and sequence 
dependentEEEd. Based upon these views the wedge model 
offered the very first explanation of bending induced by 
A-tracts by postulatina-that stacking in ApA steps is in- 
trinsically non-parallelE3. Modified versions of this the- 



ory accounted for a substantial part of available exper- 
imental data, with good scor es . p£ ^curvature prediction 
from sets of fitted wedge anglesEJii-IEil At the same time, 
clear experimental counter-examples exist where b ending 
could not result from-simple accumulation of wedgesPjEZl. 
The junction modettEl better than other theories ex- 
plained experimental data on gel retardation of curved 
DNA. It originated from an idea that a bencLshould occur 
when two different DNA forms are stackecO. If poly-dA 
double helix had a special B' form as suggested by some 
datao the helical axis should be kinked when an A-tract 
is interrupted by a random sequence. In turn, the X- 
ray data are best interpreted with an alternative theory 
that postulates that bending is intrinsic in^mos-LDNA 
sequences except A-tracts which are straightEirEjO. An- 
other interesting model attracted attention in the recent 
years, namely, bending by electrostatic forces that re- 
sult from neutralization of phosphates by .solvent cations 
trapped in the minor grooves of A-tractsa. This prob- 
lem is of general importance because the accumulated 
large volume of apparently paradoxical observations sug- 
gests that some essential features are still unknown that 
may be essential for the fine structure and the biological 
function of the DNA molecule. 

One of us has recently proposed a new hypothesis of 
the physical origin of intrinsic bends in double helical 
DNAEH According to it, the sugar-phosphate backbone 
in physiological conditions is slightly compressed, that 
is the equilibrium specific length of the corresponding 
free polymer in the same solvent is larger than that in 
the canonical B-form. Therefore, the backbone "pushes" 
stacked base pairs, forcing them to increase the helical 
twist and rise while the stacking interactions oppose this. 
As a result, the backbone increases its length by deviating 
from its regular spiral trace and wanders along the heli- 
cal surface causing quasi-sinusoidal modulations of DNA 
grooves. Concomitant base stacking perturbations result 
in macroscopic static curvature when certain properties 
of base pairs alternate along the sequence in phase with 
the helical screw. 

Drew and Travers (1984, 1985) apparently were the 
first to notice that narrowing of both DNA grooves at 
the inner edge of a bend is a necessary and sufficient 
condition of bending, and that an unusual local groove 
width should be accompanied by structural perturbations 
beyond this region. They, and later Burkhoff and Tul- 
lius (1987), considered the preference of narrow and wide 
minor groove profiles by certain sequences as the pos- 
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sible original cause of this effect. Sprous et al (1999) 
proposed a similar idea within the context of the junc- 
tion model. In a certain sense, the compressed back- 
bone theory continued the same line of thinking. Unlike 
other models, it naturally explains well-known environ- 
mental effect upon the A- tr act .-curvature , notably, its 
reduction with-tfimneraturcr^r tE3 and addition of dehy- 
drating agentsH'EH (see discussion in Mazur (2000)). 
This theory certainly needs further examination and crit- 
ical comparison with other models in both calculations 
and experiments. 

Conformational modeling earlier helped to shed light 
upon many aspects of the above problems. Construction 
of spatial DNA traces from local wedge parameters com- 
bined with Monte Carlo simulations of loop closure was 
applied to check different hypotheses and tp_ftstimate lo- 
cal bending angles from experimental datalij'fEI. Energy 
calculations revealed that bending may be easier at sQpp 
dinucleotide steps and in certain specific directionsEJo, 
with experimental sequence effects reproduced in some 
remarkable examplesO. DNA was shown to have local 
energy minima in fc$f&, conformations corresponding to 
the junction modeE3't£l. All atom Monte Carlo calcula- 
tions showed that narrowing of the A-tract minor groove 
with a few NMR-deripjd restraints may be sufficient to 
provoke the curvature£3. 

The simplest set-up for modeling DNA bending is to 
take a straight symmetrical double helix and let it bend 
spontaneously with no extra forces applied, that is due to 
generic atom-atom interactions. This "naive" approach 
has recently become possible owing to the progress in 
methodology of-molecular dynamics (MD) calculations 
of nucleic acidsES, which was demonstrated by success- 
ful simulations of several curved and straight DNA frag- 
ments in realistic .environment including explicit water 
and counterionst3'c2l. The character of bending qualita- 
tively agreed with the theories outlined above so that 
none of them could be preferred. Thorough discrimi- 
native testing would require more extensive sampling of 
bending events, which should become possible in future. 
Detailed structures otsliort A-tract fragments have also 
been studied by MDElM 

The major obstacle in free MD simulations of intrinsic 
curvature is the limited capacity of sampling of bending 
events. The physical time of transition between straight 
and bent conformations may be too long for a statistically 
significant number of such events to be accumulated in 
simulations. Moreover, experimental effects may not ap- 
pear during infinitely long MD because models are never 
perfect. To circumvent these difficulties, we employed 
a different strategy. We first looked for, and found a 
short A-tract motif that could reproducibly induce sta- 
ble bends in DNA during a few nanoseconds of MD with 
a simplified model of B-DNA. We used this motif to con- 
struct longer double helices with intrinsic curvature in 
silico and we could increase the length of DNA fragments 
in calculations to 35 bp, which makes possible a direct 
comparison with experiments in vitro. 



The two 35-mer DNA fragments we study here have 
identical base pair composition and differ only by their 
sequences. The first fragment is the designed A-tract 
repeat while the other sequence is "random" . All MD 
trajectories start from canonical straight A- and B-DNA 
conformations. For the A-tract DNA fragment they con- 
verged to a single statically bent state with planar cur- 
vature towards the narrowed minor grooves at 3' ends 
of A-tracts. The magnitude of bending is close to the 
experimental estimates. The random fragment was not 
straight as well, but its curvature was much less signif- 
icant and less planar. In gel migration assays the two 
molecules produce well- resolved distinct bands, with the 
A-tract sequence demonstrating a reduced mobility char- 
acteristic of curved DNA. These results suggest that the 
intrinsic DNA curvature reproduced in calculations cor- 
responds to the experimental phenomenon. The bending 
dynamics qualitatively agrees with the compressed back- 
bone theory, but it cannot be accounted for by other 
models. 



MATERIALS AND METHODS 
Calculations 

Molecular dynamics simulations have beetL-parformed 
by the internal coordinate method (ICMD)c3'ES-jnclud- 
ing special technique for flexible sugar ringsE3, with 
AMBER94ocZI force field and TIP3P wateiH All cal- 
culations were carried out without cut-offs and bound- 
ary conditions. The time step was 10 fseCjr-iThe so- 
called minimal model of B-DNA was usedu3'E3. It in- 
cludes only a partial hydration shell and treats counterion 
and long range solvation effects implicitly. Advantages as 
well as limitations of this approach have been reviewed 
elsewhereEa. The model has no other bias towards bent 
or non-bent conformations except the base pair sequence. 

The starting fiber A- and B-DNA models— were con- 
structed from the published atom coordinates^. The hy- 
dration protocols were same as beforepl with an identi- 
cal number of explicit water molecules in A- and B-DNA 
starts. Programs Curves,Lj XmMol,E3 and Mathematica 
by Wolfram Research Inc. were employed for graphics 
and data analysis. 

The two 35pb DNA fragments are referred to below as 
At and nAt, for the A-tract repeat and the non-A-tract 
DNA, respectively. For both fragments two long MD 
trajectories were computed starting from either A- or B- 
canonical DNA forms. These four trajectories are re- 
ferred to as At-A, At-B, nAt- A, and nAt-B, respectively, 
where the last character indicates the starting state. All 
trajectories were continued to 20 ns except At-B which 
was stopped at about 12 ns because it had clearly con- 
verged long before. 
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Experimental Sequence 



CAAAAATGTCAA^VAAATAGGCAAAAAATG . . . 3 



A 

Bend centre 
35bp A-tract repeat 



AAAATAG GCTATTTTAGGCTATTTTAGGCTATTTT- 3 ' 



3' -TTTTATCC 3ATAAAAICC 



ATAAAAICC 



Reference random fragment 

5' -TTAGATAGTATGACTATCTATGATCATGTATGATA-3' 
3' -AATCTATCATACTGATAGATACTAGTACATACTAT-5' 

FIG. 1. Construction of 35 bp double stranded DNA frag- 
ments. The top sequence with the boxed heptamer mo 
AAAATAG is taken from the trypanosome kinetoplast DNA 
The A-tracts are numbered and their centers are separated by 
approximately 10 bp. The reference random fragment has the 
same base pair content as the 35-mer repeat, but its sequence 
has been manually re-shuffled to exclude any A-tract motives. 



RESULTS AND DISCUSSION 

Construction of DNA Fragments 

Figure [t] explains how the two DNA fragments used 
in our study have been constructed. The A-tract mo- 
tif AAAATAG originally attracted our attention, in MD 

1, which 
The 35 



simulations of the natural DNA shown in Fig. 
is the first curved DNA locus studied in vitn 
bp A-tract fragment was constructed by repeating this 
motif four times and it had to be inverted to make the 
two DNA termini .symmetrical. Such inversion should 
not affect bendingpj but is essential for simulations be- 
cause the 3'- and 5'-end A-tracts may represent qualita- 
tively different boundaries. In repeated simulations with 
this and similar A-tract fragments, the static curvature 
emerged spontaneously and. it became more evident as 
the chain length increased^! To obtain a reference non- 
A-tract DNA, we have re-shuffled manually base pairs 
of the A-tract repeat. We preferred this randomized se- 
quence to commonly used GC-rich straight fragments in 
order to keep the base pair content identical and reduce 
the noise that could cause small variations in gel mobility 
and hide the subtle differences we were going to detect. 



Oligonucleotides and construction of 5'-labeled DNA 
probes 

The sequences of the 35 nt-long synthetic oligonu- 
cleotides used here are shown in Fig. [I]. The dou- 
ble stranded DNAs were obtained by annealing of the 
two complementary oligonucleotides, one of them la- 
beled with T4 polynucleotide kinase and [ 32 -P]-ATP. The 
annealing was carried out by incubating the oligonu- 
cleotides (300 nM) for 3 min at 80° C in 20 mM Tris- 
HC1 (pH 8.0), 400 mM NaCl, 0.2 mM EDTA and then 
allowing them to cool slowly. 

Gel mobility assays 

Mobility of the DNA fragments was analyzed in 16% 
gels (acrylamide to bis-acrylamide, 29:1) buffered with 90 
mM Tris-borate, 1 mM EDTA, pH 8.6. Gels were pre-run 
under constant power until stabilization of current. End- 
labeled DNA in a buffer containing 20 mM Tris-HCL, 50 
mM NaCL, 7% glycerol, pH 8.0 and bromophenol-blue 
was loaded onto the gel. The electrophoresis was per- 
formed under constant power and constant temperature 
of 8°C. The dried gels were exposed to storage phos- 
phor screens and visualized on a 400S Phosphorlmager 
(Molecular Dynamics). 



Spontaneous Development of Curvature in 
Simulations 

All four trajectories exhibited stable dynamics with 
DNA structures close to the B form. Table | shows pa- 
rameters of the final Ins-average conformations. They 
all have remarkably similar helicoidals corresponding to 
a typical B-DNA. For example, the average helical twist 
estimated from the best-fit B-DNA experimental valuesE^I 
gives 34.0 ± 0.2° and 33.8 ± 0.2° for the A-tract fragment 
and the randomized sequence, respectively. At the same 
time, the rms deviations from the canonical structures 
vary more significantly. 

As shown in Fig. g, during the first few nanoseconds, 
the rmsd from the canonical B-DNA quickly leveled at 
around 4 A in all four trajectories. For the A-DNA start 
this corresponds to a rapid transition to B-form with re- 
duction of rmsd from the initial 10.7 A. The subsequent 
dynamics is remarkably different for the At and nAt tra- 
jectories. In At-A and At-B, after some delay, the rmsd 
value drastically increased and stabilized at a higher level 
of around 6 A. The traces of the bend angle and the axis 
shortening indicate that this was a transition to a signifi- 
cantly larger curvature. In contrast, for nAt trajectories, 
Fig. H exhibits only fluctuations at roughly the same 
level as in At-A and At-B before the transition. 
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TABLE I. Some structural parameters of staadard and computed DNA conformations. The helicoidals are the sequence 
averaged values computed with program CurvesEa. All distances are in angstroms and angles in degrees. 





Xdisp 


Inclin 


Rise 


Twist 


RMSD vs A-DNA 


RMSD vs B-DNA 


A-DNA 


-5.4 


+ 19.1 


2.6 


32.7 


0.0 


10.7 


B-DNA 


-0.7 


-6.0 


3.4 


36.0 


10.7 


0.0 


At-A 


+0.1 


-4.0 


3.5 


34.2 


11.6 
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At-B 


-0.4 


-5.2 


3.5 


34.5 
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nAt-A 


-0.1 


-4.2 
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34.3 
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11.2 
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FIG. 2. The time variation of some parameters that characterize the overall DNA shape. The plates are grouped in raws 
for the same trajectory and in columns for the same parameter. The first column displays the non-hydrogen atom rmsd (in 
angstroms) from the fiber canonical B-DNAEj. The second column shows the bend angle in degrees. The last column shows 
the shortening, that is the excess length of the curved DNA axis with respect to its end-to-end distance. For example, 10% 
shortening means that the end-to-end distance is 10% shorter than the curved trace. The traces were smoothed by averaging 
with a window of 75 ps in At-B and 150 ps otherwise. 
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FIG. 3. The time evolution of the overall shape of the helical axis in At-A and nAt-A. The axis of the curved double helix is 
computed as the best fit common axis of coaxial cylindrical surfaces passing through sugar atoms, which gives solutions close 
to those produced by the Curves algorithmic. The two surface plots labeled X and Y are constructed by using projections of 
the curved axis upon the XOZ and YOZ planes, respectively, of the global Cartesian frame shown in Fig. ^a. Any time section 
of these surfaces gives the corresponding projection averaged over a time window of 400 ps. The horizontal deviation is given 
in angstroms and, for clarity, its relative scale is two times increased with respect to the true DNA length. Shown on the right 
are the corresponding views of the final Ins-average conformations. The AT base pairs are shown by thicker lines. 
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The origin of this difference is analyzed in Fig. ||. It 
displays dynamics of the overall DNA shape by using two 
orthogonal projections of the helical axis. A planar bend 
would give a plane in the Y projection and a curved sur- 
face in the X projection. A sharp increase of curvature 
in At-A after the 13th nanosecond is evident. Analogous 
event occurred in At-B after about 3 ns. In agreement 
with Fig. ||, the two nAt surfaces show fluctuations with 
amplitudes similar to those during the first 13 ns of At- 
A. This pattern probably corresponds to a generic type 
of dynamics characteristic of arbitrary 35-mer DNA frag- 
ments. 

Comparison of the three columns of plots in Fig. |^ in- 
dicates that fluctuations usually occurred simultaneously 
in all three parameters, which means that the bending 
dynamics makes a major contribution to the rmsd from 
B-DNA. Its values shown in Table | are actually much 
larger than would be for straight conformations with the 
same helical parameters. For instance, the rmsd between 
the At-A and At-B structures in Table [j] was 2.3 A only 
because, as we show below, they were bent in the same 
direction. 



Convergence of Trajectories 

A DNA molecule with detectable static curvature can 
either have a minimum of potential energy in a bent state 
or its energy valley should have a special shapp_such that 
a bent form has larger conformational entropyo. In both 
cases such state represents a free energy minimum where 
MD trajectories should be trapped. The question is, how- 
ever, how long a real MD trajectory should stay in a bent 
conformation to be representative. Some experiments 
suggest that bending dynamics in DNA fragments of only 
100 base paips.may involve relaxation times longer than a 
microsecondE3 : E£l, therefore, no practical procedure exists 
to prove rigorously that computed conformations are rep- 
resentative. Nevertheless, if several trajectories converge 
to the same state from very different starting points, one 
can argue that this state is an attractor in the confor- 
mational space, which is a necessary condition of the 
static curvature. The reciprocal convergence of trajecto- 
ries starting from canonical A- and B-DNA, therefore, is 
a very important aspect of these simulations. The mini- 
mal B-DNA model is not expected to give stable A-DNA 
structures and we did not try to equilibrate the initial 
A-DNA states. The start form the A-form is important 
because it provides an independent dynamic assay with 
a very different entry to the B-DNA family, which al- 
lows one to verify convergence of trajectories to specific 
conformations. We analyze separately two levels of struc- 
tural convergence. 
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FIG. 4. (a) Geometric constructions used for evaluating 
the DNA bending. The two coordinate frames shown are the 
global Cartesian coordinates (OXYZ), and the local frame 
constructed in the middle point of the curved rBNA axis ac- 
cording to the Cambridge convention (O'JKL)EJ. The curve 
is rotated with two its ends fixed at the Z-axis to put the mid- 
dle point in plane XOZ. The bending direction is measured by 
angle ip between this plane and vector J of the local frame. By 
definition, this vector points to the maj|«j DNA groove along 
the short axis of the reference base paiiO. Consequently, the 
zero ip value corresponds to the overall bend towards the mi- 
nor groove in the middle of the DNA fragment, (b) The time 
evolution of the bending direction as measured by the ip angle 
in plate (a) (in degrees). The traces have been smoothed by 
averaging with a window of 75 ps in At-B and 150 ps other- 
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Overall shapes 

The rmsd comparison between At-A and At-B is shown 
in Fig. S. It clearly demonstrates that At-A and At-B 
trajectories managed to come very close to each other 
even though their starting points were significantly sepa- 
rated in conformational space. The initial rmsd of 10.7 A 
between the canonical 35-mer A- and B-DNA forms even- 
tually went down to as low as 1.3 A. The final fall of the 
rmsd occurred when the curvature drastically increased 
(compare Figs. ^ and ||). Moreover, during the last 
nanoseconds the bending direction was virtually iden- 
tical in At-A and At-B and essentially fixed at around 
90° (see Fig. |]b), which explains the origin of the black 
rectangle in the upper right corner of Fig. |^. This di- 
rection corresponds to bending towards the minor groove 
at approximately three base pair steps from the middle 
GC pair (see Fig. [|a), that is at the 3' end of the third 
A-tract in Fig. [l]. 

The nAt trajectories exhibited qualitatively different 
features. The rmsd comparison of any two long intervals 
of nAt-A and nAt-B gives fluctuations between 3 and 6 
A without any clear time trend. Figure |^ shows that 
the rmsd from B-DNA also fluctuated between 3 and 6 
A and that it correlated with bending parameters. As 
seen in Fig. || the molecule really was not straight. Ac- 
cording to Fig. the bending directions in nAt-A and 
nAt-B were well defined but slightly different. They nei- 
ther diverged nor converged, remaining at around 100° 
from each other. The molecule shows no signs of slow 
straightening which would give a decrease of fluctua- 
tions in Fig. and an increase in scattering of directions 
in Fig. |Jb. All this suggests that bent shapes are fa- 
vored over straight ones, and that there are many stable 
bends, with transitions between them being too rare to 
be sampled by our simulations. 

Groove profiles and local structures 

Dynamics of the minor groove profiles is shown in Fig. 
^. There are evident qualitative resemblance as well as 
some subtle differences between these four surfaces. In 
At-B, the characteristic regular groove shape has estab- 
lished early, with significant widenings in the three zones 
between the A-tracts. The very left widening is some- 
what different probably because it occurs between anti- 
parallel A-tracts. In At-A, the profile strongly changed 
at the beginning, but also established by the end of the 
10th nanosecond. Although the final At-A and At-B pro- 
files are not identical, they are clearly similar, with good 
correspondence of local widenings and narrowings. 

The two nAt surfaces show little similarity between 
each other, but qualitatively their shapes are not very 
different from those for the A-tract fragment, with modu- 
lations of similar wavelengths and amplitudes. This looks 
somewhat counter-intuitive because, in experiments, reg- 



ular oscillations of the pinor groove widths are observed 
only in A-tract repeatsEj, and this structural periodicity 
is certainly related to that of the sequence. However, such 
behavior is exactly what one should expect if the waving 
of the backbone results from its intrinsic compression. In 
this case, the groove modulations should occur regardless 
of the base pair sequence and their characteristic wave 
lengths should be determined by the backbone stiffness 
as well as overall helical pitch and diameter. This ex- 
plains why the waves in the left-hand and the right-hand 
plates in Fig. ^| have roughly similar scales, even though 
only the A-tract sequence is periodical. In experiment, 
however, such modulations can be observed only if their 
phases are fixed in time, which is the case of A-tract re- 
peats. For random sequences, like the one we use as a 
reference, the fine structure should be smoothed out on 
averaging over the whole ensemble. 

Figure |^ compares Bi/Bn backbone dynamics in the 
two At trajectories. There are many similarities in dy- 
namics as well as in the final configurations. The conver- 
gence is better near both ends and within A-tracts. The 
dissimilar distributions of the conformers in the middle 
corresponds to the difference in minor groove profiles in 
Fig. |^. In A-tracts, the Bn conformers are very rare 
in T-strands and tend to alternate with Bi in A-strands. 
Figure || compares local helical parameters in the last av- 
erage structures. Only the Buckle and Propeller traces 
exhibit large scale modulations phased with the helical 
screw. All parameters strongly fluctuate and these fluc- 
tuations are apparently chaotic with rather dissimilar 
phases in the two structures. 

Coupling between the levels 

Figures || and ^ demonstrate that At-A and At-B ar- 
rived at the same statically bent state. This dynamics 
contrasts those of the two nAt trajectories and it strongly 
suggest that the curved DNA shape of the A-tract frag- 
ment is an attractor of trajectories with a meta-basin 
of attraction comprising both canonical A and B DNA 
forms. Figures [||8| show that the bending convergence is 
accompanied by some clear trends in local conformational 
dynamics. These local features are probably coupled to 
bending, however, a close look reveals that this coupling 
is very loose. The convergence of the minor groove pro- 
files in Fig. H is at best qualitative. Figure fj] indicates 
that active backbone dynamics continued after the cur- 
vature has established and that one can pick up rather 
different distributions of conformers from the ensemble 
of bent structures. The noisy traces in Fig. || obtained 
by averaging over two similarly bent ensembles suggest 
that the helical parameters are far from being constant. 
The natural conclusion follows that convergence of the 
bending dynamics does not require unique specific local 
conformations, i. e. that the bent state is microhetero- 
geneous. 
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At-B 

FIG. 5. A 2D density plot of the rms difference between At-A and At-B. Conformations spaced by 2.5 ps intervals were first 
averaged over of 50 and 25 ps intervals in At-A and At-B, respectively, and the resulting structures compared between the 
trajectories. Darker shading implies smaller rmsd values. The lower left corner corresponds to the initial structures, that is the 
canonical A- and B-forms, with rmsd about 10.7 A. The shaded rectangle in the upper right corner demonstrates convergence 
of the two trajectories to the same bent state. The values above 4 A are not shaded whereas in the darkest zones it falls down 
to 1.3 A. The black vertical band at approximately 10 ns indicates that At-A shortly visited the final state 3 ns before the 
definite transition. 




FIG. 6. The time evolution of the profile of the minor groove in the four trajectories. The surface plots are formed by 
time-averaged successive minor groove profiles, with that on-the front face corresponding to the final DNA conformation. The 
groove width is evaluated by using space traces of C5' atomsS. Its value is given in angstroms and the corresponding canonical 
B-DNA level of 7.7 A is marked by the thin straight lines on the faces of the box. The sequences are shown for the corresponding 
top strands in Fig. |l| with the 5'-ends on the left. The A-tracts are underlined. Note that the groove width can be measured 
only starting from the third base pair from both termini. 
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FIG. 7. (a) Dynamics of Bi and Bn backbone conformers in At-A and At-B. The Bi and Bn conformations are distinguished 
by the values of two consecutive backbone torsions, e and £. In a transition, they change concertedly from (t,g~) to (g~,t). 
The difference ^ — e is, therefore, positive in Bi state and negative in Bn, and it is used as a monitoring indicator, with 
the corresponding gray scale levels shown on the right. Each base pair step is characterized by a column consisting of two 
sub-columns, with the left sub-columns referring to the sequence written at the top in 5'-3' direction from left to right. The 
right sub-columns refer to the complementary sequence shown at the bottom, (b) Comparison of the final distributions of 
Bi and Bn backbone conformers in At-A and At-B shown in the same way as in plate (a). 
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The Magnitude and The Character of Bending in 
the A-tract Repeat 



The experimental magnitude of bending caused by A- 
tracts was earlijesjestiroated by several groups with differ- 
ent approachesEirEill£3'L3. The reported bend angles were 
between 11° and 28° per A-tract, and 18° .is presently 
considered as the most reasonable estimated. The cur- 
vature somewhat varies with the base pair sequence and 
depends upon environmental conditions such as the tem- 
perature, the concentration of counterions etc. Although 
in calculations all these details cannot yet be properly 
taken into account a quantitative comparison with ex- 
periment is instructive. 

When the curvature has established, that is after 13 
ns of dynamics in At- A and after 3 ns in At-B, the bend 
angle oscillated around 60° (see Fig. ||). In the consec- 
utive Ins-averaged conformations its value was between 
42° and 74°, with the average of 54° for 16 such struc- 
tures. This value corresponds to 54/4=13.5° per A-tract, 
that is close to the lower experimental estimate. A larger 
value of 54/3=18° results, however, if one assumes, as 
suggested by some experimental observationsE]l£3 that 
the A-tracts are straight, and that the bending actually 
occurs in the three zones between them. Yet another 
estimate is obtained from the increase of bending with 
respect to the shorter 25-mer fragment studied earlierEl 
It appears that one additional A-tract and junction zone 
increase the overall bend by 20-22°. We see that the 
magnitude of bending in simulations is rather close to 
experimental estimates, and that the agreement is better 
if the curvature is really localized in the junction zones 
between A-tracts. 

Figure || presents a closer look at how the local cur- 
vature is distributed in the last 1-ns average structure of 
At-B. The total bending angle is about 50°. Three zones 
contribute more than other to the overall bend. The 
two junctions between A-tracts 2, 3 and 4 are bent in 
an identical direction which is close to that of the whole 
structure. Together they contribute around 40° to the to- 
tal bend, which is the largest local positive contribution. 
In contrast, the strongly curved fourth A-tract makes a 
negative contribution because its direction diverges by 
more than 90°. The third A-tract is virtually straight. 
Finally, A-tracts 1 and 2 and the junction zone between 
them exhibit a smooth curvature with a stable "good" 
direction and contribute the remaining 20° of the total 
bend. 



AAAATAQGCTA' 




TTTTATCCGMAAAATCCX3ATAAAATCCX3ATAAAA 



FIG. 8. Sequence variations of helicoidal parameters in the 
last Ins-average structures of At-A and At-B. The sequence 
of the first strand is shown on the top in 5' - 3' direction. 
The complementary sequence of the second strand is written 
on the bottom in the opposite direction. All parameters were 
evaluated with the Curves programtj and are given in degrees 
and angstroms. At-A - solid line, At-B - dashed line. 
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FIG. 9. (a) The last 1-ns average structure of At-B shown in the XOZ projection according Fig. ^a. The AT base pairs 
are highlighted, (b) The quantified distribution of curvature in the structure shown in plate (a). The local bending angle 
is evaluated by moving a sliding window along the helical axis. The window size was 3 base pair steps, with the measured 
values assigned to its center. The sequences of strands are given as in Fig. [j] with the A-tracts underlined and numbered, (c) 
Dynamics of local bending in At-B. The surface plot is formed by time-averaged successive profiles like that in plate (b), with 
the front face of the box corresponding to the end of the trajectory. 
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The foregoing analysis certainly is not free from pit- 
falls. For instance, the apparent smooth curvature can 
result from time averaging of several alternative local 
bends. Nevertheless, Fig. g indicates that there are zones 
in this DNA fragment that are bent more than other and 
that two such zones are distinguishable between A-tracts. 
Figure ^|c displays the local bending dynamics in At-B. It 
is seen that the main features noticed in plates (a) and (b) 
were quite visible during the whole trajectory. Moreover, 
the zone between the first two A-tracts also sometimes 
carried an increased curvature. However, it would be in- 
correct to conclude that A-tracts are straight. They just 
exhibit generally smaller and distributed curvature than 
the junction zones. This curvature is usually directed to- 
wards the minor groove, therefore, it does not cancel out 
in averaged structures. 

The foregoing pattern agrees qualitatively with the re- 
cent NMREa and X-ray dataE3 as well as-the. character of 
bending earlier observed in calculationsEj'El Many ear- 
lier reported X-ray structures of A-tracts suggested that 
they produce an-intrinsically straight DNA compared to 
other sequenced^. Our calculations do not contradict 
these observations because the crystal A-tract structures 
should be additional^pteaightened due to special crys- 
tallization conditionscS3E3, and because a single short 
A-tract may in fact be somewhat less curved than that 
inserted in a long DNA fragment. 



Verification of Curvature by Gel Electrophoresis 

The sequence induced static DNA curvature was first 
noticed owing to reduced migration rate of curved DNA 
fragments in gel electrophoresis!!]. Later gel migration 
studies provided a, wealth of information on curvature 
in A-tract repeatso. The difference in gel mobility be- 
tween straight and curved DNA rapidly grows with chain 
length, therefore, the curvature was usually studied in 
rather long DNA fragments. Data for sequences shorter 
than 50 bp are rareEH, and, to our best knowledge, it has 
never been shown that curved and straight 35-mers could 
be distinguished. Nevertheless, subtle sequence effects in 
double stranded oligomers of .around 10 bp were detected 
with higher gel concentration!!], and one could hope that 
this would work for somewhat longer sequences as well. 
Alternatively, the effect of bending could be enhanced 
by inserting constructed fragments into a long stretch of 
straight DNA0, but this would complicate further analy- 
sis because the DNA molecules used in experiments and 
calculations could no longer be identical. 

Figure [l^ shows comparison of the acrylamide gel mo- 
bility of these two fragments. As expected, the A-tract 
repeat exhibits a reduced rate of migration. The differ- 
ence is quite significant so that the two molecules are well 
resolved both in separate lanes and when mixed in the 
same sample. Owing to the identical base pair content, 
the minor factors such as the number of tightly bound 



Gel start 




nAt At At+nAt 

FIG. 10. Gel mobility assay. The two 32 P-labeled 35 bp 
DNA constructs (At and nAt) were electrophoresed in 16% 
plyacrylamide gel buffered with Tris-borate, pH 8.6. The gel 
was dried and autoradiographed. The lanes labeled At, nAt, 
and At+nAt correspond to the A-tract repeat, the random 
sequence, and their mixture, respectively. Bands assigned to 
each DNA fragment are marked by arrows. 
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covmterions and water molecules is reduced here to the 
possible minimum and, most probably, the observed dif- 
ference is entirely due to the curvature in the A-tract 
fragment. 



DISCUSSION 



tions and we believe this approach presents considerable 
interest for future studies. Sequence effects in DNA frag- 
ments of 35-50 base pairs can be probed in both MD 
simulations and gel electrophoresis. Such experiments 
are rapid and inexpensive, which progressively becomes 
the case for MD simulations as well. 



Comparison with Earlier Studies 



Comparison with Theories of DNA Bending 



To our knowledge, the only earlier successful unbiased 
simulations aimed at reproducing A-tract induced cur- 
vature_jp_DNA have been reported by David Beveridge 
groupLJO. These simulations were carried out in full 
water environment with explicit counterions. The char- 
acter of the phased A-tract bendins-appeared oscillatory 
with a period of at least 3 to 4 nscJ. Because duration 
of trajectories was only 5 ns, it was difficult to confirm 
the static character of bending and distinguish between 
essential and occasional observations. Therefore, con- 
clusions concerning applicability of different models were 
not restrictive and left room for many theories. Our sim- 
ulations have the same goal and a similar set up, but we 
use a simpler model system. The primarily long term 
interest in B-DNA models with implicit or semi-implicit 
representation of environment is connected with approx- 
imate simulations of very long DNA moleculesES. As 
shown here the minimal model can also capture, at least 
qualitatively, important sequence effects like the A-tract 
induced curvature. 

Several features in our calculations correspond well 
to those observed earlier, notably, spontaneous develop- 
ment of quasi-sinusoidal minor groove profiles in both 
A-tract and non-A-tract sequences and strong bends in 
junction zones between A-tracts. In contrast to earlier 
simulations, however, the curvature here emerged after 
several nanoseconds of dynamics and the difference be- 
tween the A-tract and non-A-tract structures did not 
reduce with time. One should note also that the A- 
tract structures computed with the minimal model are 
very close to experimental data as |*epards the helical 
pitch and the absolute groove sizesE3E3. In standard 
AMBER and CHARMM simulations, B-DNA always ap- 
pears somewhat underwound and the narrowest A-tract 
minor gropves, jcemain 1-2 A wider than in experimental 
structures2J'c2rc3. The origin of this subtle bias remains 
unclear, and attempts to reduce it have been made ip_i 
very recent modifications of the AMBER force fieldcJ 
In the minimal model, this bias was compensated by fit- 
ting reduced phosphate charges, which apparently im- 
proved the A-tract structures and stabilized the curva- 
ture. 

A significant enforcement oi-.tlie present results com- 
pared to our previous reportsEdCj consists in the direct 
comparison of bending in silico and in vitro, which be- 
came possible owing to increased length of the DNA frag- 
ment. Such a possibility is rather unique for MD simula- 



The origin of intrinsic curvature in DNA remains un- 
clear. Theories that explain it always assume some spe- 
cific balance of interactions in the DNA structure, and 
that is why these theories are perhaps more important 
than the particular role of A-tracts. The list of available 
interactions is well-known, but the question is which of 
them is the driving force. Below we briefly analyze our 
results in contexts of some theories. 



Base Pair Stacking Models 

According to any mechanism that starts from hase.aair 
stacking, like the wedge or the junction modelsc3'E£ra, a 
curved DNA must be built out of asymmetric blocks, 
with their structures determined by base pair sequence. 
The bending, therefore, must be accompanied by rep- 
etition of local structures in identical sequence frag- 
ments. This fundamental theoretical prediction fails for 
the static t^pds observed here, which confirms earlier 
conclusions£3a. The structures of sequence repeats in 
the bent state are microscopically heterogeneous and con- 
vergence to specific local conformations is not necessary 
for bending. As shown above, the A-tract trajectories ar- 
rive at a single bent state, but the minor groove profiles 
in Fig. H are only similar, not identical as well as local 
helical parameters and backbone conformations in Figs. 



and 



Counterion Electrostatic Models 



An alternative model that recently attracted much 
attention considers solvent cations trapped in-A-tract 
minor grooves as the initial cause of bendingcJ. The 
role of cousteijkwis in this phenomenon is rather 
controversial 6 iri'LJ , and a few general comments are nec- 
essary before considering our results. Because straight 
DNA structures correspond to symmetric minima of elec- 
trostatic energy bends can result from symmetry break- 
ing in the charge distribution, namely, if positive external 
charges accpntJat|C ftt-iQne DNA side it should bend to- 
wards theralll'OOO'Ej. However, the same situation is 
well interpreted by other models of bending. Namely, 
in a curved double helix, the phosphate groups at the 
inner edge must approach, which creates regions of low 
potential that should be populated by counterions if they 
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are available^. In the first case the counterion-DNA in- 
teractions are sequence specific and they cause bending. 
In the second case they are structure specific and they 
stabilize pre-existing curvature. 

Two physically different models of sequence-specific 
counterion involvement can be distinguished. In the first 
one the counterions act locally. When a counterion is 
placed in one of the DNA grooves between two phos- 
phate groups their electrostatic interaction becomes at- 
tractive, which narrows the grooveU. As in some ear- 
lier modelalj'Ej, the global curvature results from a gen- 
eral mechanical link between groove deformations and 
bending. In contrast, the second model is purely electro- 
static. HereJ-he minor grooves of A-tracts act as flexible 
ionophorescJ'Lj and trap counterions. Since in phased 
sequences they occur at the same DNA side the double 
helix bends towards them to relax the long range phos- 
phate repulsion at the opposite side. 

The second model employs the general idea initially 
proposed for protein DNA ipteractionsO and confirmed 
experimentally for free DNAO. However, it qualitatively 
disagrees with a cornerstone experimental observation 
concerning the A-tract induced bending, namely, that 
an A-tract can be characterized by a definite bend an- 
gle regardless of its length and the distance from other 
A-tracts. When the length of an A-tract exceeds one heli- 
cal turn both sides of the double helix appear neutralized. 
As a result, the curvature should decrease in the series 
(Ai 2 Ng) n — (Ai4N7) n — (Ai 6 N 5 ) n because the length of 
the non-neutralized N-tracts is reduced, and furthermore, 
in sequence (Ai 6 N 5 ) n the bend angle per A-tract should 
be drastically reduced with respect to that in (A 6 N 5 ) n , 
for example, because the distance between the repulsive 
N-tracts is increased. These predictions apparently dis- 
agree with the experimental trendsEZl although additional 
experiments are perhaps necessary to check them. 

The first model cannot explain the origin of the A- 
tract curvature because-only multivalent counterions can 
cause significant bendalj whereas bending is commonly 
observed in buffers containing EDTA and other chelat- 
ing agents. Also, the optimal counterion position for this 
type of bend is at the entrance of the groove and not 
inside, therefore, it cannot be both strong and sequence 
specific. The last argument agrees with the recent MD 
studies of correlations between the minor groove width 
and positioning of counterions. Notably, there is no such 
correlation when only counterions interacting with bases 
are consideredEx In contrast, a correlation, exists for 
counterion positions at the groove entrancdia. The last 
observation corresponds to the structure specific bind- 
ing better than to the sequence specific one. Structure- 
specific interactions can explain all experimental results 
concepaiiig-the preferential binding of counterions in A- 
tractslIaaH, which makes such data intrinsically neutral 
as regards different models of bending. 

In our calculations all counterion effects are considered 
non-specific, and the results obtained indicate that mod- 
ulations of DNA grooves and static bending are phys- 



ically possible without breaking the charge symmetry 
around DNA. Although simulations alone cannot prove 
the real mechanism all experimental and computational 
observations taken together suggest that solvent counte- 
rions are hardly responsible for the intrinsic curvature in 
DNA, which by no means questions their important role 
in DNA structure and function. 



Compressed Backbone Theory 



The main idea of the compressed backbone theor 




was outlined in Introduction. All seemingly paradoxical 
MD observations from which it originally emerged are 
confirmed here. Notably, this theory predicts that, with 
any base pair sequence, the backbone stiffness should 
cause smooth modulations of DNA grooves. The helical 
symmetry becomes broken with the base pair stacking 
perturbed, which creates regions of intrinsic curvature. 
In a "random" DNA, the local curvature changes its di- 
rection with time because groove widenings and narrow- 
ings migrate slowly along the double helix. As a result, 
the generic DNA appears straight on average although 
it is curved locally. In sequences where certain base pair 
properties strongly alternate, the phases of backbone os- 
cillations appear fixed. In this case the local curvature 
can sum up to give macroscopic static bends, as in A- 
tract repeats. This theory considers a macroscopically 
curved DNA as an "idioform" characterized by topolog- 
ical attributes, rather than a structure with fixed atom 
positions. These attributes are the bend direction and 
the phase of groove modulations. The microheterogene- 
ity of the bent state should be expected because the same 
waving backbone profile is compatible with many alter- 
native local conformations. 

A competition between the stacking interactions and 
the backbone compression postulated by this theory-is 
characteristic of physical systems called frustrated!^. 
Consider the common textbook example of three anti- 
ferromagnetic spins in a triangle configuration. The op- 
timal orientation of each pair is anti-parallel, but all three 
pairs cannot be anti-parallel in a triangle. There is always 
at least one parallel pair and the ground state appears de- 
generate. Now consider a circular duplex DNA with a ho- 
mopolymer sequence. The compressed backbone causes 
groove modulations, but there are no preferable regions 
for narrowings and widenings and the ground state ap- 
pears strongly degenerate. The similarity between these 
two examples is evident. One usual physical consequence 
of frustration is very important for biology, namely, the 
possibility of a glassy state where microscopic transitions 
are dramatically slowed down. Transitions between wavy 
backbone configurations in a long DNA can be very slow 
because many groove narrowings and widenings must be 
moved concertedly. This may explain observations of 
supra-microsecond relaxation times ia_ bending dynam- 
ics of relatively short DNA fragmentsEj'El 
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Finally, the compressed backbone theory offers a new 
view of some environmental effects upon the curvature. 
Common physical factors like the temperature, counte- 
rions, and various dehydrating agents aa a lon g-known to 
change slightly the helical pitch of DNAeSHM^ which can 
reasonably be attributed to the dependence of the state 
of the DNA backbone upon the solvent screening of phos- 
phates. These same factors modulate significantly the se- 
quence specificity of nucleases probably by changing the 
shape of DNA groovesEj, a»d-p»jidjpe|--complex effects 
upon the intrinsic curvaturel^'EJ'ESEJOS. It seems wise 
to postpone any detailed interpretation of these facts for 
future studies, but one can just note that with intrinsic 
frustration outlined above a very small change in the par- 
tial specific backbone length can induce significant global 
changes in the DNA structure. 



Possibilities of Experimental Verification of 
Backbone Compression 

The the compressed backbone theory does not give 
simple rules for a priori calculation of curvature in any 
sequence. Nevertheless, it offers some predictions that 
can be checked in experiments. It suggests, for example, 
that the A-tract curvature can be relaxed by introducing 
single-strand breaks. It seems interesting also to examine 
the possible relationship between the backbone compres- 
sion and supcrcoiling. There is a consensus that intj-ia- 
sic bends affect the shape of the superhelical DNAEjEZI. 
Unlike other models, however, the compressed backbone 
theory predicts that the intrinsic curvature should vary 
under superhelical stress in a rather special way. Namely, 
with a positive density, the backbone is stretched and 
the curvature of an internal A-tract repeat should be re- 
duced. Conversely, the curvature should increase when 
the superhelical density is negative. Diekmann and Wang 
earlier observed that ti^e A-tract structure changes un- 
der superhelical stress^!, and their approach may serve 
for a more specific experimental verification of the above 
predictions. 
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